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SUMMARY 


: ye 

Α method is presented for determining the time-dependent flow over 
a rectangular wing moving with a supersonic forward speed and undergoing 
small vertical distortions expressible &s polynomials involving spanwise 
and chordwisé distances. The solution for the velocity potential is 
presented in & form analogous to that for steady supersonic flow having 
the familiar "reflected area” concept discovered by Evvard. Particular 
attention is paid to indicial-type motions and results are expressed in 
terms of generalized indicial forces. Numerical results for Mach numbers 
equal to 1.1 and 1.2 are given for polynomials of the first and fifth 
degree in the chordwise and spanwise directions, respectively, on a wing 
having an aspect ratio of ἃ, 


INTRODUCTION 


One of the basic problems arising in the analysis of wing flutter 
boundaries is the calculation of the aerodynamic forces on wings under- 
going small but arbitrary spanwise and chordwise distortions. When the 
wing aspect ratio is large (actually, when the distance between spanwise 
nodal lines is large), these forces are usually estimated by some strip 
theory in which the loading on each spanwise section is approximated 
from that on a two-dimensional wing having the same chordwise distortion. 
This report is concerned with low-aspect-ratio rectangular wings for 
which tip effects are important and the full three-dimensional theory 
must be used. 


The exact linearized solution for the forces on thin rectangular 
wings (limited, however, to the range where effective aspect ratio 


(/ M2-1 A) is > 1) traveling at supersonic speeds has been presented by 
both Gardner (ref. 1) and Miles (refs. 2 and 3) in terms of multiple 
integrals involving arbitrary surface undulations. However, the use of 
such solutions in evaluating, numerically say, the forces induced by 
specific wing distortions still presents some difficulties. It is the 
purpose of this report to discuss certain techniques that can simplify 
the labor involved in these calculations and to present numerical tables 
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for the forces induced by a class of surface deformations, a class gen- 
eral enough to represent the first few mode shapes of rectangular 
plates. 


Mathematically the problem is to find and analyze a solution to the 
four-dimensional wave equation 


Pax + Pyy + Paz n Puro: = 0 (15) 


(where a, is the speed of sound, t' is the time, and x,y,z are space 


coordinates) that satisfies the appropriate boundary conditions. ‘The 
particular form of the solution to be analyzed differs from those pre- 
sented by Gardner and Miles but its development is based on the method 
due to Gardner. 


Hadamard (ref. 4) studied a generalized form of equation (la) in 
which the number of dimensions was arbitrary. His solutions to these 
generalized equations are fundamentally different, depending on whether 
the total number of dimensions is odd or even. In fact, the methods 
Hadamard developed apply directly only to equations for which the total 
number of dimensions is odd. Solutions for the even cases (such as 
eq. (la)) are determined by a "method of descent"; that is, the solution 
for the next higher odd-dimensioned equation is found and then reduced 
by (made independent of) one dimension. It is apparent, however, that 
such & technique is in itself by no means unique. Thus, Hadamard found 
the solution to equation (la) by descending from a solution to the 
equation 


l a 
Ye + Pyy + Paz * θεὲ -E Peres = 0 (1b ) 


but there are many other partial differential equations and groups of 
partial differential equations governing a five-dimensional (x,y,z,t,t) 
space all of which satisfy equation (la) in a plane € = constant. 
Gardner discovered a set of equations containing equation (la) in a 

t£ = constant plane which are simpler than equation (la) in that solu- 
tions could be found and adapted to the boundary conditions for time- 
dependent motion by methods well known to aerodynamicists who have 
studied the flow about wings in steady supersonic flight. This is the 
essential part of Gardner's contribution and it represents the tech- 
nique upon which the development of the solution presented in this 
report is based. Actually, Gardner first applied a Lorentz transforma- 
tion to equation (la) and then used his method outlined above. The 
application of such a transformation is unnecessary and has the dis- 
advantage that the resulting coordinates have lost their direct physical 
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significance. We will apply Gardner's method of descent directly to 
equation (la) and then proceed to analyze the solutions so obtained. 


In order to simplify the analysis as much as possible, we will 
limit solutions to the plane of the wing, and, further, consider only 
indicial-type boundary conditions; in other words, unsteady motions in 
which the wing attains instantaneously, at time zero, a certain span- 
wise and chordwise distortion which is thereafter fixed. It is well 
known that the transient responses to these indicial motions can be used, 
in a superposition integral, to obtain responses to many other types of 
unsteady motion; in particular, responses to the harmonic oscillations 
of nonrigid wings. 


Finally, the principal interpretation of the results will be made 
in terms of generalized forces, since these can be used directly in 
either flutter or gust studies, and it will be shown that the amount of 
labor required to calculate such forces is reduced by using reciprocity 
relations derived from the general theorems presented in reference 5. 


LIST OF IMPORTANT SYMBOLS 


A aspect ratio 

Bo speed of sound 

aln amplitude of indicial-downwash distribution (See eq. (2a).) 
B(p,q) beta function (See eq. (Bl5a).) 


B,-y2(p,g) incomplete beta function (See eq. (B15b).) 


CGe ya) influence function for effect of side edge (See eq. (A10).) 


Cr, lift coefficient, lift 
doS 
Siy indicial lift coefficient due to angle-of-attack change, 
OCT, 
without pitching, UE = — 
da | a0 
Cha. indicial lift coefficient due to pitching for & wing rotating 
ac 
about its leading edge, στα Ξ 2 
δα la 
Cm pitching-moment coefficient, positive when trailing edge 


tends to sink relative to leading edge, E 
qs Sc 
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indicial pitching-moment coefficient due to angle-of-attack 
change (without pitching) measured about the leading edge, 


ο0πι 


indicial pitching-moment coefficient due to pitching measured 
about the leading edge for a wing rotating about its lead- 


OCm 


wing chord 


generalized indicial force coefficient (See eq. (36).) 


generalized indicial force coefficient (See eq. (37).) 


distance of wing camber line from z = O plane 
Mach number 


loading coefficient (pressure on the lower surface minus 
pressure on the upper surface divided by free-stream 


dynamic pressure) 


n! 


binominal coefficient, 6 = ni (nai 
m. n-m). 


dimensionless rate of pitching, 2 
O 


free-stream dynamic pressure, =PoUo? 


generalized coordinate 


generallzed force corresponding to the generalized coor- 
dinate qy 


real part of 


af (x-x, )* + (y-y)* 
N (x-x,)^ + (y*y4)* 


af (x-x, )* E BS (y-y,)* 
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X, (η) 


wing semispan 

wing area 

area of acoustic plan form 

area of reflected acoustic plan form 
agt! 


time 


wing kinetic energy : 


wing potential energy 


forward speed of wing 


(= 
92 /z=0 


vertical velocity 


Cartesian coordinates, fixed relative to the fluid at 
infinity 


coordinates with origin on center of wing leading edge 
(See sketch (1).) 


coordinates with origin on center of wing leading edge at 
time zero (See sketch (n).) 


Ω |. 


Mx + t 


B 
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α angle of attack (angle between flight path and plane of 
wing), radians 

p M^ - 1 

0 wing angle of pitch relative to horizontal, positive when 
trailing edge lies below leading edge, radians 

£ coordinate measuring fifth dimension 

Po free-stream density 

p velocity potential 

φί3) portion of velocity potential induced by sources in acoustic 
plan form 

"o portion of velocity potential induced by presence of side 
edge 

Y potential function in five-dimensional space. 

iN 
Subscripts 

A,B,C regions in an x,É plane (See sketch (d).) 

u upper side of wing, z = O+ 

1 singularity (e.g., source) position 


I,II,...VIII regions on wing shown in figure 1 


STATEMENT OF THE PROBLEM 


The Governing Equation 


Assuming a wing's vertical motion is of such a nature that the 
velocities induced in the fluid are small relative to the magnitude of 
the wing's steady forward motion, the normalized form of equation (la) 


xx yy zz jet rue (10) 
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where t = agt!, can be used as the governing partial differential equa- 
tion of the flow field. This equation applies to the determination of 
the velocity potential when the body or wing in question moves through 
the fluid, the axes remaining fixed with respect to the still fluid 
infinitely distant from the origin. For convenience we place the wing 
leading edge on the y axis at t = O and the side edge on the x 

axis. The wing flies at a constant forward (in the negative x direc- 
tion) speed so at subsequent times the leading edge lies &long the line 
x = -Mt, where M is the Mach number, and the side edge moves along 

the x axis as shown in sketch (a). 


xz - Mf 


The Boundary Conditions 77 
The fluid velocity normal to the Yj Yj 
surface of a solid moving in a friction- ΖΑ 


less fluid must be zero. If the equa- 
tion of the solid's surface is repre- 
sented by 


G(x,y,z,t!) = 0 


this boundary condition can be expressed Sketch (a) 
mathematically, in terms of the coordinate system used in equation (le), 
as 


dG, OG, ME, ὁφ OG _ 


== + = + --- --- - 


ðt! Ox Ox dy ὃν Oz Oz 


Consider a thin surface near the z = O plane. The equation of the 
camber line of this surf&ce can then be expressed in the form 


G(x,y,z,t!)- z - h(x,y,t?) = 0 


and, assuming that thickness and lifting effects can be separated lin- 
early, the boundary condition for the camber line becomes 

Oh , 9P oh , OP dh 99. 

ðt! Ox dx dy Oy dz 


If the derivatives of h with respect to each of the coordinates are 
small, the two middle terms can be neglected and the expression for the 
boundary condition reduces to 

Eh dd ei. = πε κει 
Ot Oz =, 
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We wish to simulate a rectangular wing deformed indicially by bend- 
ing in the n chordwise directions. For this purpose, on the 
portion of the = O plane occupied by the wing plan form, the vertical 
velocity, which ανν πάς the wing shape ἜΠΗ to the previous 
equation, is &ssumed to have the form 


O » στο 
Wy = l n 
x+Mt 
L 2 (E) (2) 
1 n 


where c is chord length, azn is a constant and wt and n are 
integers > 0. 


The expression (x + Mt)! is used so that for 1 > O the tangent 
to the wing camber line at the leading edge is tangent to the flight-path 
angle of the leading edge. Consider, for example, the case 2=1, n=0. 
The downwash | 


a 
Wu = -Ξ- (κ + Mt) 


represents an infinite class of surface shapes having the form 


h(x,y,t) - [(x + Mt)? + f(x,y)] (e) 


CUO 


where f(x,y) is an arbitrary function and h is, by definition, the 
distance of the wing's camber line from the z = O plane. Since, within 
the accuracy of linearized theory, the solution for the flow about the 
wing depends only upon the value of w,(x,y,t), the loading on all the 
wings represented by the above equation is the same. 


Let us inspect the two special cases 


^ 


(1) f (x,y) 


(ii) f(x,y) 


li 
O 


For case (i) 


M 
19. (2xt + Mt?) 
2οῦο 


h(x,y, 
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and the wing is a flat plate pitching at a uniform rate about its lead- 
ing edge which is following the flight path 


(h) = _ Bight t? 


LE 2cUg 


as shown” in sketch (ο). Hence, αὖ time t the tangent to the flight 
path of the leading edge is 


, le 2 " 


ΠΠ... Dim E gistribution 


flight path of ~ S 
leading edge . 
ag ο) 
Ν Γ σ 
Sketch (ο) b 


à(h)../dt' azot 


E ς 
Uo 


The slope of the leading edge of the plate at the same time is 


à) Mor 
Ox LE c 


and the two slopes are seen to be equivalent. 


For case (ii) 


a 
h(x,y,t) - 2 19. (x + Mt)” 


cU, 


‘The z scale in both sketches (b) and (ο) is purposely distorted 
in order to make the drawings clear. A basic assumption used in setting 
up the boundary-value problem, by means of which the loading was deter- 


mined, was that the surface of the wing must remain near the z-O plane. 
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and the wing is a plate which obtained a sudden parabolic camber at 
t = 0, a shape it maintained thereafter as shown? in sketch (c). 


σε ΞΟ at leading edge 


Sketch (c) 


The problem is linear, so it will be sufficient to determine & 
solution for arbitrary l and n, and then add results for any combina- B 
tion of terms as desired. Thus, the complete boundary conditions to 


be studied are 
t n 
CE) (3) . 
Z=0 


over the wing plan form, and, since the loading is zero over the remain- 
ing portion of the plane 


Wu (x,y,6) = a 


99 
ot 


= 0 off the wing (2b) 
z=0 


since the loading is given by 


A E 
do | UgM «Ot /z=0+ 
SOLUTION FOR THE POTENTIAL 


Figure 1 shows the wing plan form on the surface of which the 
potential is required, together with the system of axes; &lso, traces 
in the z = O plane of the wave system set up by the indicial motion of 
the wing are indicated. The wave pattern for only two edges is shown; 


“See footnote 1 on p. 9. 
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the flight speed is supersonic so the trailing edge has no effect on 
the velocities induced over the wing surface, and the results are valid 
(in their entirety) only for fA » 1, so the opposite edge either has 
no effect or one that can be incorporated by simple superposition. 


The wave traces divide the wing area into several regions, indicated 
by the Roman numerals, in each of which the analytical formulation for 
the potential is different. Region I consists of that part of the wing 
where the effect of neither the side edge nor leading edge has yet been 
felt. In region II, the side-edge influence is acting (the line y= t 
is the trace of the starting cylindrical wave from the side edge y = 0) 
but not the leading edge. Region III is the part within the starting 
cylindrical wave from the leading edge, but outside the influence of the 
side edge. This region, and region V, are further subdivided for reasons 
that will appear later. Region IV is a compound region; potential there 
can be found by adding the potentials for regions ΤΙ and III and sub- 
tracting the potential for region I. Region V consists of the portion 
of the wing within the spherical wave originating at the wing corner. 

The flow over the part of the wing comprising regions VI and VII has 
reached & steady state relative to & point on the wing, and the poten- 
tial there is just that for the corresponding parts of & rectangular 
ving with the proper downwash distribution in steady motion. Finally, 
region VIII is again a composite region, its potential being the sum of 
potentials for regions III and VII less the potential for region VI. 


All the regions just listed, with the exception of region V, are 
actually governed by the three- (total) dimensional wave equation and 
the potential therein could be obtained by methods applicable to this 
simpler equation. However, in this report we shall present a unified 
approach and the problem will be solved by the same method in all 
regions. 


Review of Kirchhoff's Formula 


The solutions developed in the subsequent sections are more clearly 
interpretable if they are compared with certain known results that have 
already been determined for the indicial motion of nonlifting wings with 
symmetrical thickness distributions or lifting surfaces with all super- 
sonic edges. The purpose of this section is simply to review briefly 
some of these latter results. 


As in steady-state wing theory, there is a formula for time- ” 
dependent flows that relates the velocity potential to & distribution 
of time-dependent sources and doublets over & certain region in the 
wing plane. This formula is due to Kirchhoff, and some of its aero- 
dynamic uses are discussed in reference 6. Kirchhoff's result is 
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immediately applicable in the study of unsteady lifting-surface problems 
when the potential can be represented by sources alone, that is, when 
the upper and lower surfaces of the wing do not interact, as is the case 
in regions I, IIT, and VI of figure 1. 


Kirchhoff's formula for source distributions can be written 


Q(x,y,0,t) = = 57 ή el dx, dyi (3) 
Da, 


where 


το = (x - xi) + (y - σι) 
The brackets on w indicate that the retarded value is to be taken 
[Wu] = Wax, ,Y 1» t-ro) 


and Sa indicates that the region of integration is the acoustic plan 
form corresponding to the event (x,y,0,t). These concepts are discussed 
at length in reference 6. 


As has been pointed out, equation (3) holds for each of the regions 
I, III, and VI, but the area of integration Sa differs considerably 
from one of these regions to another. Consider, for example, the deter- 
mination of 1 for region III, denoted PITT: Part of the boundary of 
the acoustic plan form Sa is found by eliminating T between the 
equation of the leading edge, x, = -MT, and the expression 


(x - x)? + (y - y,)* = (t - m* 


which gives the outer boundary, at "time" t, of all the disturbances 
that, operating at "time" T, can produce an effect at the point (x,y). 
This boundary is the ellipse 


(£ x, - m) t(y-y)?5t (ha) 


where 


. 
*- A. ve meam ASA P P er ni E em ag 
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If the point (x,y) lies within the cylindrical wave from the leading 
edge, that is, -t< x < t, the ellipse of equation (4a) comprises only 
part of the acoustic plan form, the remainder being bounded by so much 
of tbe circle 


(x -x)^-(y-y^ = tf (4b) 


as lies on the wing at time zero. Sketch (d) shows the three possible 
acoustic plan forms for points in region III. The limits for the three 


types are 


(1) t>x>0 
(11) 0>x>-t/M 
(111) -t/M » x > -t 


(ral t- = 


(ex J} (y- y) 5 1* 


Sketch (d) 


and these correspond to the subregions Illa, IIIb, and IIIc identified 
in figure 1. Using equation (3), we can write the potential in, say, 
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region Illa as 


y+t x+a/ t^-(y-y,)? 
"ta " 7 Bx / ayy "n El dx, + 


yt ti? Xi(y-yi) ἘΝ 
= n dy; a dx, (5) 
y-4 t?-x? x- t*-(y-yi)* 


where 


Xi(y = yi) = : | κα E tm - TA 


Gardner's Method of Descent 


Equation (1) governs a four-dimensional x,y,z,t space. Our object, 
of course, is to find for this equation a solution that satisfies the 
boundary conditions in the z = O plane as specified in equations (2a) 
and (2b). Obviously, we can always construct a space of more dimensions 
governed in an arbitrary way except that it must satisfy equation (1) in 
an X,y;Z;t hyperplane. Then, if a solution in this higher dimensional 
space which satisfies equations (2a) and (2b) in the x,y,z,t plane can 
be found, it represents for £ (the additional dimension) equal to some 
constant the solution to our problem. This characterizes the method of 
descent. It is not obvious, of course, that such a method leads to any 
simplification; but, with & proper choice of the governing equation for 
the new space, such a possibility always exists. 


There are examples where various applications of this method have 
proved to be useful. Hadamard's use of the method, mentioned in the 
introduction, is classical. A simple application of his method is the 
derivation of the velocity potential for & source in & two-dimensional 
supersonic flow field. This potential field (which amounts to a step 
function, the step occurring at the Mach wave) is easy to derive if one 
considers a three-dimensional field with a line of sources normal to the 
free stream and uniform in strength. The two-dimensional field mentioned 
above follows immediately by descent. 


In other examples the additional dimension is measured with imagi- 
nary numbers and the additional law for the extended space is the 
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requirement that the functional dependence on the resulting complex 
variable shall be analytic. The method of descending in the latter 
cage is associated with the study of analytic continuation. In partic- 
ular, Riesz!s method (discussed in ref. 7) for solving equation (1) 
illustrates these concepts. 


Gardner's method for solving equation (1) is to define a five- 
dimensional space in which a potential function Y is governed by the 
. equations 


Ver = Yyy = Lm = 0 (6b) 


and show that solutions to equations (6) in this space are general 
enough to contain general solutions to equation (1) in a plane 

E = constant. We shall, therefore, proceed by analyzing these equa- 
tions and eventually let € approach & plane in vhich the boundary 
conditions of equations (2a) and (2b) are satisfied. For convenience, 
the latter plane is taken to be the & = O plane. 


Since equations (6a) and (6b) are linear, a number of possibili- 
ties exist for the choice of the dependent variable ¥(x,y,z,0,t). 
Aside from the more obvious choice Y(x,y,z,0,t) = Q(x,y,z,t), where 9 
is the velocity potential of equation (1); for example, one could let 
V(x,y,z,0,t) = Qu (x,y,z, t) or again, Ve (x,y,2,0,t) = p(x,y,z,6). These 
various choices amount only to relatively minor differences in the 
detailed technique of the subsequent analysis. If, in imposing the 
boundary conditions of equations (2), one is to use only source-type 
solutions for both equations (68) and (6b), the last choice is suffi- 
cient. Therefore, set 


| 2 V(x,y,z, t,t) | = (X,Y, 2» Ὁ) (T) 
οξ E=0 


Now differentiate equation (6a) with respect to z and set? z=0. 
“Tt can be shown that the solution satisfies the equation 


SS [RS [nemen] = SS tenen 
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Defining 


W(&,x,y,t) x et 5 (8) 
z= 


equation (6a) can be expressed in the form 
Wit - Woe - Weg 50 (9) 


and the boundary conditions for equation (9) are given directly by equa- 
tions (2). Thus on the wing 


2 n 
9ξίε-ο — Ozlg.o 5 j 
and off the wing 
A D (x,y ,0,%) = 0 (10b) 
dt leno 


Assuming equation (9) to have been solved for the boundary condi- 
tions given by equations (10), we return to the second of the set of 
partial differential equations (6), specifically, 


Vee - Yyy - Ven = 0 


From equation (8), it is seen that the solution to equation (9) 
yields the result 


oy 


= known function of y,& on the wing 
Oz 70 


Further, the boundary conditions for the original problem in (x,y,z,€,t) 
space require that  ὍὌθ an odd function with respect to z, and con- 
tinuous across the z = O plane except over the wing plan form. Thus € 


3M 
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must be zero for z = 0 except over the wing plan form. The continuation 
of this condition into (x,y,z,&,t) space then implies, according to equa- 
tion (7), that off the wing 


S oa 
abl s . 


Hence, both the second p&rtial differential equation and its boundary 
conditions are identical in form to the first set given by equations (9) 
and (10), respectively. Applying equation (T) to their dual solution, 
we obtain the desired result 


E Z] .- Q(x,y,0,t) 
dE E =O 


for the potential on a rectangular wing (with A >1) in supersonic 
unsteady motion. 


The General Expression for the Potential 


* 


The method outlined in the preceding section will now be applied 
to obtain integral expressions for the potential in any region of the 
rectangular wing shown in figure l. Consider first equation (9) for 
W(&,x,t). This equation is the same partial differential equation as 
that which governs supersonic steady flow. Further, the boundary values 
in the €&,x,t space are identical to 
those representing a thin planar wing 
in a steady supersonic flow. Since the 
Mach number in the steady-flow analog 
is N 8, the equivalent plan form of 
this wing (shown in sketch (e)) is & 
sweptforward wing tip having all super- 
sonic edges (1.e., the component of the 
free-stream velocity normal to all edges 
is supersonic). 


Since all edges of the equivalent 
wing plan form are supersonic, the 
solution for W can be written imme- 
diately in terms of "sources" only, 
their strength being given by equa- 
tion (10a). Thus, by analogy with | Sketch (e) 
the well-known results of supersonic 
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wing theory, we have 


Wy, (Xx7+M6, , y) dx, dt, l 
553b; = - TI 
" Ñ A eroe (t-t,)*- ES- (x-x 1) i 


where T is the area on the wing cut out by the forecone from the 
point (£,x,t), see sketch (e). The analytic form of W will differ 
considerably in each of the three regions above me equivalent wing 
shown in sketch (f). 


The value of W given by 
equation (11) now becomes a bound- 
ary condition for the solution of 
equation (6b). Thus, over the 
portion of the z=0 plane for which 
y ΣΟ, E20, the variation of 


Ν 2 is now known and for 
2170 
i, y <0, &>0O the condition 
e 5 = 0 applies. (These condi- 
á 9ξ z=0 


tions are still not sufficient to 
determine a unique solution unless 
the further restriction is imposed. 
that the loading falls to zero as 
the edge y Ξ O is approached, 
i.e., as y—>0+.) Again we observe 
that these boundary conditions and 
the partial differential equa- 
tion (6b) are identical to those 
studied in connection with a sta- 
tionary planar wing in a supersonic 
stream. As shown in sketch (f), solutions from the +t,x,& space above 
the £t = O plane are referred to as Wa, Wp, and Wo, depending on the 
relation between x and ἔ ina t = constant plane. Sketch (g) shows 
the five different boundary-value problems formed by the various com- 
binations of Wa, Wa, and Wo occurring along constant x lines in the 
x,& plane; and the corresponding regions in figure 1 for which each 
applies. Each of these five problems is directly analogous to the 
boundary-value problem encountered in steady-state lifting-surface 
theory, of & planar, rectangular lifting surf&ce in & steady supersonic 
stream. The "leading edges" of these analogous rectangular plan forms 
lie along the lines δ. = t, ZW = A t^-x* or E, = tm, depending on the 
value of x, and tbe "side edge" lies along the ins y ^0. Hence, by 
means of this steady-flow analog, we can immediately write the solution 


NN 


SS 
ES 


Sketch (f) 
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to equation (6b) in the form 


, Y Y(x,y,0,E,t) = 
y ? W(x,y,,5,,t)d5,dy, 
Yj (12) 
7, AN TSS οἱ 


J| (&£-54)* - (y-y,)* 


where only the area of integration σ 
must be discussed. 


J, 
(i) x2; for regions LT, 


Two possibilities exist for the 
shape of o. First, if the point 
E, y lies to the right of the dashed 
lines in sketch (g), which in the 
‘analogous steady-flow problem repre- 
sent the traces of the Mach cones 
from the leading-edge tips, σ is the 


triangular area shown (for region 
(11) Osx=t; I, III&) in sketch (h) part (i). If 
é 


however, ἔ, y lies between this line 
and the side edge, y = 0, o is the 
trapezoidal area shown (for region Va) 
in sketch (h) part (ii). The latter 
is & well-known result used in steady 
supersonic lifting-surface theory and 


S 
(iii) -t/Msxs0; I,E, 


(i) O<x<t , y»! 


N 
(yv) -Mtsx<-?; LT, WT 


(ii) O«x«t , y<VPFx* 


Sketch (g) Sketch (h) 
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first developed by Evvard (ref. 8). The division of the five kinds of 
problems illustrated in sketch (g) into the final twelve, represented by 
the regions in figure 1, is brought about by the various combinations of 
Wa,» Wa, and Wa that can occur in the area o as the point €, y assumes 
all necessary values on the wing. 


When y has been determined, the potential in the physical plane is 
found by equation (T), or, combining equations (11) απᾶ (12), 


Q(x,y,0,t) = 


1 lim 3 d$ ay Wy (xy Mt, , y, )üx, dt, 
Ἐν. τς 
AA a ex) 


A detailed analysis of equation (13) for a point x,y,t in region 
Va οἳ figure 1 is given in Appendix A, and & study of this analysis 
enables one to write the results for all regions without difficulty. 


Interpretation of the Results 


The results of the rather involved analysis given in Appendix A 
can be interpreted in terms of the known solutions for simpler boundary 
conditions. These latter solutions have already been reviewed in a 
previous section in which it was shown that the potential on a lifting 
surface with all supersonic edges can be written in the form 


Wy ldx 
Q(x,y,0,t) = - + n 7 [Wy Jax, dy; 
2x To —— 
Sa 


From Appendix A it is found that the potential at a point on a rectangu- 
lar lifting surface can always be expressed as the sum of two parts 


p(x,y,0,t) = φί3) (x,y,0,t) τ p (2) (x,y,0,t) (114) 
where l 
(1) 2 1 [vu Jdx, dy; 
o" yt) = - d E T DRM (158) 
Sa 
and 


p (x,y,0,) =- z JI C(x, ,y,)dx,dy, (15b) 
ο | 


ς 
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The value of C(x,,y,) is given by equation (Α10) in Appendix A and the 
areas of integration, S4 and Soe, are illustrated for the various regions 
I through VIII in figure 2. 


Let us first inspect equations (15) in light of their possible 
analogy with the familiar solution for the steady-state, rectangular 
lifting surface. If & rectangular wing having arbitrary twist and camber 
is placed in a steady supersonic flow, the solution for the potential on 
its surface can also be expressed as the sum of two parts 


e(x,y,0) = e (x,y,0) - φί3) (x,y,0) (16) 
where, if 
το = (x - x,)%- BÜ(y - σι) 
sy, 2-2 ff mA - (178) 
S 
and s 
leo) =~ S (170) 


These equations can be construed in 
the following simple way: Equation 
(17a) represents the potential induced 
at x,y,0 by & distribution of 
sources over the wing plan form, each 
source having a strength proportional 
to the local streamwise slope of the 
upper surface. The area 394, as 
shown in sketch (i), is the portion 
of tbe wing within the Mach forecone 
from x,y,0. Equation (17ο) has a 
similar interpretation; it also rep- 
resents & distribution of sources 
over the wing, each having & strength 
proportional to the local slope of 
the upper surface. But the area of 
integration Sə is now that portion 
of the wing within the Mach forecone 
from the point x,-y,0; that is, 
within the cone which forms a mirror 
image of the physical Mach forecone 
in the vertical plane containing the 


dec o — η 


———— ---- A OS 
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wing's side edge. The potential pl2) (x 7,0) represents the difference 
between the potentials for a wing with a vertically symmetrical thick- 
ness distribution and & surf&ce with no thickness having the same shape 
as the upper surface of the nonlifting wing. 


Let us return now to equations (15). Just as in the steady-state 
case, p1) (x,y,0,t) represents the potential induced at x,y,0 bya 
distribution of sources (see eq. (3)) over the wing plan form, each 
proportional to the local slope of the wing, but now, since the wing is 
in motion, with the added condition that they be local slopes at the 


P7 (ys) d 


Sketch (3) Sketch (k) | 


appropriate time. The area Sa, shown in sketch (3), is just the 
acoustic plan form defined earlier in the discussion of equations (3) 
and (4). Physically, Sa represents those points on the wing from 
which disturbances can, at the time t, influence the flow at x,y,0. 
It is the generalization, in the stationary coordinate system, of the 
wing &re& bounded by the Mach forecone. 


The relation between oJ (x y ,0,t) and (2) (xy ,0, €) is simllar 


to that between their steady-state analogs. Thus, again, o (2) (x ,y,0,t) 
represents the difference between the potentials for an uncambered non- 
lifting wing and & lifting surface having the same shape &s the top of 
the nonlifting wing. A more striking similarity lies in the relation 
between Sg and Se. 


We have already seen that Sa is the acoustic plan form, and, as 


ed 
el Lcd 
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(In other words, Sa is the acoustic plan form for the event X,y,0,t, 
and Se is the acoustic plan form for the event x,-y,0,t.) Physi- 
cally, Se represents the portion of the wing's lower surface contain- 
ing disturbances which can, at the time t, influence the flow at 

X,y,0 on the wing's upper surface. At this point the similarity 
between the steady and unsteady solutions ends since the influence of 
the slopes in the reflected plan form is not the same as it is for the 
Slopes in the basic acoustic plan form; the influence in the former 

case now being given by the integral C(x,,y,) defined in equation (A10). 


One ean shov, by simply referring the results given in equations 
(15) to a coordinate system fixed on the wing, that equations (15a) and 
(159) are identical, respectively, to equations (178) and (17b) when 
they apply to regions VII «and VI in figure 1; regions in which, for 
indicial-type motions, the flow is steady relative to the wing. Hence, 
equations (15a) and (15b) extend Evvard!s "reflected area" concept to 
&ll parts of & rectangular wing in supersonic unsteady motion.* 


THE GENERALIZED FORCES 
Review of Lagrange's Equations of Motion 
In order to define more clearly the subsequent concepts and nota- 


tion, we will briefly review Lagrange's equations of motion as applied 
to distorting wings and will examine a simple application to a rectan- 


gular wing. 


Lagrange!s equations are usually written 


zd. 9D. Οὔ. ΟΣ Gu Be να. α (18) 


dt' 0g, da, 0d. 


where 

T kinetic énergy of the wing 

U εναν energy of wing 

Qr 8 generalized (external) force 


dy  generalized coordinate 


“It is of further interest to notice that equation (15b) can be 
reduced to a double integral involving wul bsy) by using, for example, 
the transformations € = χι Mt, and T-i-t, and integrating in the τ 
plane. 


—— Mo —ÀÁ pa — M —— M — — — — —7 +a — 
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In the present application q, is the amplitude at a given time of a 
polynomial measuring h, the. vertical displacement of the wing's camber 
line from the z = O plane. Thus, relative to an Xs, ya coordinate 
system that is fixed on the wing, see sketch (1) 


h(x yyt") = * q (8! )Ρε Gr σα) (19) 


The wing's kinetic energy can be written 


T = ff 5 h°n(xg,Vg) ἀκοᾶγα (20) 
S 
where m 15 the wing mass per unit plan-form area. Using equation (19), 
we find 
d oT » 
PX dg J | Ῥε(αο ισα) Ρο Cra sva nC γα) ataye 
dt 98. 
S S 
(21) 
o 
αν, 


The potential energy is usually difficult to evaluate analytically. 
However, it can often be determined 
experimentally (as will be seen) by 
measuring the frequeneies of the free 
vibration modes.. For the present 
assume that the wing is & homogene- 
ous plate of constant thickness. 

The potential energy for such & 

wing can be expressed as (ref. 9) 


U-2 y {Faem B ES 


ces Ó } ἀχράνο (22) 


Sketch (1) 


4M 
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which leads to the equation 


O-P,, OPP 
v?p ep. = δα pu r 8 
E Dw ff WW ; dI 27^ 


O?P, 95Ρ IP OFP 
ο ος Εις τς Β ) |esete (23) 
E a Oe  Oxg9yg Xas 


where u is Poisson's ratio, V? = ὀξ/οχα + d /ογα”, and 


2(Young's modulus)(plate thickness ) 


3(1 - pî) 


D = 


Now, if the generalized coordinates have been normalized so that 
each measures the amplitude of a free vibration mode, all terms in equa- 
tions (21) and (23) involving the integral of the product of Py and Pg 
are zero. Assuming, henceforth, such normalization, we can write 


S 2 5.. 12 O^P, OP, 
dy [ τ Ῥγ-(κα»Υα)π(κα»γα)ἀκαθγη + Da, / E {ο Pp) -2(1 -u) | z = 
S [ο OX Oyg 


BP. Ne 
( : ; }ρόκοάσο- Qr; r=1,2, +... (24) 
OXaOy s 


Finally, dividing through by the coefficient of d, and expressing a 
generalized force as the integral over the wing plan form of the product 
of the rth mode shape and the loadings® y (^p)g induced on the wing by 
each of the mode shapes considered, we find 


A 
do È fp Prora (E) dx4dys 
9 g Ο 5 


d, + quu, A --------ς -- --- (25) 
i d μ ῬοΞ(Χα»Υω)π(Χα»ΥΩ) ἀκαᾶγα 


where wy is the frequency of the rth’ free vibration mode. 


We will write (Ap) qo (Ap/a,) where q, is the free-stream 
dynamic pressure. This 0 posible without a confusion of notation since 
the generalized coordinates are expressed as 4124244». + - . and 
exclude the term qo- 
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if the free-mode frequencies are experimentally determined, equa- 
tions - such as equation (23) - giving the wing's potential energy, 
never have to be evaluated. Further, in such cases, equation (25) 
applies to quite general wing structures with varying density. Usually 
in the application of equation (25), one uses the actual frequency ὦ. 
of the free mode but, in evaluating the aerodynamic forces, uses an 
analytical expression that only approximates the rth mode shape. Let 
us examine the generalized force term in equation (25), taking, for 
simplicity, only one term of the sum; 


ar = ao ff 7G) ς bis (26) 
o 


According to what has gone before, the mode shape polynomial P,(xs,ys) 


has the form 
Bre) -Y (Y (27) 


while (Ap/a,). is the loading coefficient corresponding to an indicial 
dePlection (see previous section on boundary conditions) 


a = a (ay (9) (E) | (28) 


which gives a vertical velocity distribution 


va = Yo ας) (28) (2) (29) 


Now a generalized inđicial force coefficient can be defined as follows: 


Re 30 SO wa] en 
S 


(The calculation of these quantities fe) will be elaborated in the 


next section.) Since the generalized force Q, is intended to apply 
to any motion, not necessarily indicial, it is necessary to apply 


Duhamel's integral to the indicial force coefficient e(t); thus, 


in 
t' fi (1?) 
Sr = doS er tglt! = τη) qs (1) d (31) 


——————————— 
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ΑΒ an example, consider now the simple one degree of freedom 
vibrating plate illustrated in sketch (m). The plate is fixed to the 


[ 


Sketch (m) 
wall and restrained along its leading edge. The mode shape is assumed 


to have the form 
2 2 
x 
h = q,(t,') (2) (38) (32) 


so for a plate with uniform density and thickness 


O ς ^ " 4 
m 7 Ya f dxa Pr^(xa,ya) = eS (8) 
Ο 


-8 


Equation (25) now becomes 
ie ' 
.. 2 S 
τω = — 


For this case, we have the generalized indicial force coefficient 


1911 
f a(t ), and so 


t! £19(q1) 
Qa, = qetse) le [ q (t-T') E oan dT! (34) 
dt! 
A q, (1) 


Therefore, equation (33) can be written 


. 25 y tU rom) 
ito a (2) TA a (61-11) | 22— | art (35) 
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The Generalized Indicial Force Coefficient 


It is clear from the previous section that a study of the dynamic 
behavior of rectangular wings moving at supersonic speeds can be carried 
out if one can obtain values of the generalized force coefficient, 


ey), as defined by equation (30). We will now show how these values 
can be obtained from the solution to the aerodynamic boundary-value 
problem represented by equation (14). 


Το was convenient in developing 
equation (14) to use a coordinate 


Xy = Xe ή system - x,y,z,t - which was fixed 

Js 7 Yq in space so that the left edge of the 
237 24, wing moved along the x axis as 

ta = ty shown in sketch (a). On the other 


hand, in studying the dynamic problem 

15 was more convenient to use an 
1,20 XarY a 410, system which was fixed 
in space so that the wing's spanwise 
center line moved along the x, axis, 
see sketch (n). Let us first con- 
sider the problem of transferring the 
results in terms of the x,y,z,t 
coordinates to the x,,y,,Z2,,t, 
system. 


The indicial force coefficient 
Patt) is defined as follows: 


Sketch (n) 


; c-Mt 8 3 g 1n 

Dry 1. x+Mt y Ap 

Pagl*") "το x f ay ( c ) (5) (2) vn 
-Mt O 


In order to transfer the axes from the set shown in sketch (a) to the 
more convenient set of sketch (n), so that mode shapes are symmetric or 
asymmetric about the wing's spanwise center line and the force coeffi- 


cients denoted p. can be determined, we »roceed as follows. First, 


the loading coefficient for & wing in the (x,y) system with downwash 
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given by 


«σος GE) er) of OOO 


is obtained. This loading coefficient can be written as a sum: 


(E Tl y (a) E (0) 


μ-ο 


Now the quantity a is defined as 
c-MG 
i 
in 1 E ολη X (ΔΡΧ 
f. = — dx, dy, RS = 

JE  2sc sep ο ο do 

- -S 

c-Mt 


28 d g in 
ES a ay( 2) (22) (22) 
2se ς ς do 
~Mt O 


This last integral can be written as 

c-Mt 5 
στ fT QI) Ge" 

JE  2sc ο ο do 
-Mt O . 
gin | g ν n E 
ώση creer ἢ er y" à 
2 ν 2 p 2 se 
V=0 μ-ο 


ο-Μο 


f = fey Θ᾽ 


-Mt 


- — ——— ma 
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By using equation (36) we find 


abet y arg Σο» ΘΩ͂ rm α 


where all forces are responses to a unit indicial disturbance. Note 
that if equation (37) is applied in the case of a wing cantilevered on 
a wall, both n and g must be even in order to satisfy the boundary 
conditions of reflection in the wall. 


By superimposing boundary conditions &nd their resulting solutions, 
one can further show that the value of " given by equation (3T) is 
valid for all reduced aspect ratios PA greater than 1 in spite of the 
fact that the value of P given by equation (36), as it stands, 
applies only to wings for which BPA is greater than 2. 


Given get), one can determine the generalized force associated 


jg 
with the generalized coordinate q, by means of the superposition inte- 
gral as illustrated by equation (34). 


Details of Calculation 


The details of actually evaluating the indicial force coefficients 
from the solution for the potential presented in the first part of this 
report are discussed in Appendix B. Considerable labor is involved in 
such calculations, and an attempt was made to discover recursion formulas 
by means of which certain derivatives, for the rectangular wing, could be 
expressed as combinations of others. This attempt was successful and 
yielded the following results 


Consider equation (36). Integrate the x integral in this equation 
by parts, setting 


in 
8 
u(x) v ys £P ay; dv(x) = (x+ Mit)? dx 
S do 
Then, since by equation (BT) in Appendix B 


στ -- ---- ————— --- ----.-- 
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one finds 
in i l-1,n o 
= ——P - F. a 
F^ πμ ο 8 jeg (384) 
Inspection of equation (37) shows that the same relation holds for the 
generalized indicial force coefficients "s that is, 
in i l-1,0 _ ¿l-1,n 38b 
de^ jel {πο g iu em 
From this relation, it is seen that only the forces P need be 


determined by integration; the forces for higher values of the index 1 
can be found by combination of results for different values of the mode 
shape index  j. 


As a simple illustration of the results presented so far, we can 
calculate the indicial force derivative for the cases l=n=g¢=0, 
j20,1. The case ἡ = O corresponds to the indicial lift coefficient 
for a flat, sinking, rectangular wing, and the case for j = 1 corre- 
Sponds to the indicial pitching-moment coefficient for the same wing. 
Since n= g = 0, equation (37) gives 


Thus, with ἡ = O and identifying -200/Uo as angle of attack a, one 
finds from Appendix B 


1 oo ἃ to Mto 1 
στι = -—— f -π|ὶ--π(ι--ο 0 < to < — 
ii 200/00 99 M | 7 ( e ΠΌΤ 


= 4 ji ο. o + : cos^l(M- β: 5ο) + Y to* - (1- ΚΗ : 


O 


1 1 2 1 1 
alm οί > |} = Σ 


safi mn to > — 
— M-1 


Next, with j = 1, and using Cn! to designate the pitching moment 


TET ae — A ty Án mae A n má E MÀ R 9 at üt n [9 — — — d — 
— — À— -— ~e 
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measured &bout the leading edge of the wing, 


1 oo 2 NT 2 
ong! = -(- x 5e 4( 5 to”) E (Mé+ 1)to J} 


2 
2 {2 | € = 2) cos” Oe 4 5 cos”1 (M- B%,) + 


LMt 
? fe? - Mo) | - = ÉS 3 to - (M-1) tos |) 


These expressions agree with those given by Miles in reference 2. 


The above results can be used to demonstrate the usefulness of 
equation (38a). Taking j=n=g=0, 2 = l in that equation gives 


10 ΟΟ OO 
Foo oo ~ Fio 
or, for the present case, 
10 00 00 
oo 7 foo D fio 


which represents the equality 


q Cr, + Ong' 


that is, the lift coefficient for & pitching wing equals the sum of the 
lift and pitching-moment coefficients of a sinking wing (primes indicate 
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the wing is pitching about and moments are measured about the wing lead- 
ing edge). Hence, 


aie 1 1 2 1 Mty-1 M -1 4. RE 
- 21 2[ (+B to?) 00s e το (M- Bto) + 


to 
3 - Mto 2| 1 L 
to^ - (1 - Mt ) FEE. 3to - 3(M-1)t eu 
2 O 6A Mtl O O 

-ote |} -ᾱ- < to < -+ 
M+1 ~ Mel 

-24 ου to > 1 

β 3βΑ — M-1 


A further application of equation (38a) provides the pitching- 
moment coefficient for & pitehing flat rectangular wing. Thus, with 
l=j=1l,n=8=0, equation (38a) gives 


which becomes 


and so 


-> —— aera, 
ree -- -- -------ϕ- ---- mm ee rr ee ag A gg r ~ 
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From equation (B21) in Appendix B it is found that 


£99 sd 

20 5ο kh fy 8 to 5 

200 80 M 3 7 Mto ) - Fo.) a «3x! OXto X v 

Uo Uo 
_ 4 ffilm? y Mto-l 1 M anaal 
"Sell poop pt 
l+Mto+H (M2) t0 fo .. 2 1 
M d tof - (1 - Μες) ] Aue - 
9 j il OhA| Mel j 
8 4 1 L 
M - A ee 
"Viu tas LL 
813 ` IBA κ Mel 
Combining, we find 
1. 2 2+Mto το 8 - Mt + M(M? + 3)t,° O t -L 
Cmg M 3 19A O O < O < M+1 
"TS 2 uà 2+ Mt? -1 Mto - 1 ΟΜ ~] 
= pa | cos Eni xh (M - Bĉto) + 
8-Mt ¿-(M*+2)t,* Pn uec 
ο ες E G, -- | 3. + 8t, - 6(M-1)t,? 
9 o -( to)” | πιστα, ο (M-l1)to" + 
M-1)? κο Hl. 45 2 

as) |} ΜΙ 


Another relation among the generalized indicial forces "E can be 


derived by means of the reciprocity relations given in reference >. The 
details of the derivation are given in Appendix C and there results 


νων (1) A ) (1) (1 = (39) 


μ-ο H=0 
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Equation (39) can be used in two ways; one, as a means for checking 
the internal consistency of a set of calculated generalized indicial 
forces, and the other, as a means for expressing a given force in terms 
of & set of others. 


Consider, 86 an example of the former use, the case for which 


on og 
og on 


From equation (37) bis can express this relation in terms of the calcu- 
lated quantities Po thus 


Poi eno as 
yay (s Ye aX wm 


u=0 
If now n = 1l, g = 3 the P relation results 


CEOE CEDROS 
ΟΣ 


which provides & useful check on the computed quantities. 


o 
ro 
O FP 
Na 
+ 
PES 


Next let us solve equation (39) for a given force. Perform the 


sum operation 
- J 
je) 
zo Γ 


on both sides of equation (39), and reverse the order of summation on 
the left side. There results 


Por feq RR e 


HO j=} 


—— —— —Ó — 
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The inner sum on the left can be evaluated. Thus one has 


xP = [1 - (1-x)]? = TS (2) (ix), 


μ-ο 


| 
SOS 
Fore or)( 


rzo 


Equating coefficients of x, 


LAO ae 


and equation (40) becomes 
i d j/3 t 1,38 
n , QU 
Pjg 7 | ( -1) 53 (-1) (1 Pn (41) 
j=0 ΗΞΞΟ 


CONCLUDING REMARKS 


A method is presented for evaluating the generalized forces on & 
rectangular wing flying at supersonic speeds and having an aspect ratio 
such that ΒΑ 21. The generalized coordinates used to define the wing's 
behavior are the amplitudes of downwash distributions expressed in terms 
of polynomials in x and y, the chordwise and spanwise directions, 
respectively. 


Numerical results are presented in tabie I for generalized indicial 
forces on a wing having an aspect ratio of l and flying at a Mach number 
equal to 1.1 and 1.2; the polynomial coverage being O € 1 € 1 and 
O <n $5, where w~ xy’. 


Ames Aeronautical Laboratory 
National Advisory Committee for Aeronautics 
Moffett Field, Calif., June 30, 1954 
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APPENDIX A 
EXPRESSIONS FOR THE POTENTIAL 


In order to write the expressions for the potential in all regions 
shown in figure 1, it is sufficient to derive in detail only that for 
region V. Having carried out this analysis, one can determine the expres- 
sions for potential in other regions without difficulty. 


Consider, therefore, equation (13) and let o and T apply to region 
Va. First, it is necessary to determine the potentials Wa, and WR in 
the t,x,t space. From equation (11), in conjunction with sketch (f), 
it is found that 


2 Eo tal (x-x4)^* E? 
xt E, (x-x,) "+ By Wy (x1 + Mti ,yi)dt; 


E E [ dx, ut αρα σα ό 
i VOEE 
x-/t^-E, O (t-t,) ~ 6) > - (x-x1) 
(A1) 
ο t-,/ (x-x 7E = 
Wye -2 [ ax, / ΟΡ 
π 
χη (£1) x, /M y (t-t,)- - E c -(x-xi)* 
/ 2 ἔ 2 . / τ 2 2 
TS P Wy (Xy + Mb yy) at, (A2) 


d. 
ef ar, ER ET 
ο ο (t-t,) - & ? - (x-x4) 


where 


X,(§,) = iG Ε Jin? 597) 


With the values of W given in equations (Al) and (A2) it is pos- 
sible now to solve equation (6b) for Y, sketch (g) giving the required 


μα μου AAA ο A AAN 
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data in the £,y plane. Thus, if ΒΞ = (E-t,)? - (y-y,)? 


y t ~E+y+t t W 
W A 
¥(E,x,y,t) = - i ay, f ag, 2-5 / ay, | db, = - 
E+y-t 8a (y-y,) y E-(y-y1) 
2 V 
Wa-W 
τ dy, dt, = E 


, Ery- A t^-x? E+(y-y,) 


A γρ 


E Wp-MA 1 f^" 
ah af a] 


P t-(y-y1) E+y-t 
t -E+t-y t 
Wa 1 WA 
ds, + ij dy, dé, E + 
E(y-y,) 2 E+(yty,) 
O t?-x? 
ijf ay, as, Hawa, 
Exy- t^-x* E+(y-y,) 
2. 
ρα VP ty 
z | : ES. (A3) 


ο E+(y+y,) 
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Now apply the operation of equation (7) and the potential Qy 15 


given by = 
y t Es ΥΠΟ t E 
TE ^ A Α 
ντα fa af o uc. 
y-t Υ-γι Y -(y-y,) 

i Wi atu). AAA ror. λα 
1 ay, a BAS, " τ £ MAA. 
y -/ ta y=Y1 e 7 (y-y,) Ri 

ο pex 
EW - 
H ay, $ aby LT Ji y, F d$, Zw: Ya) 
y-t  Yy-Yi; + y-tt-x* τσι E 
t-y WA | να t 
BY. 
J dyi εντ... Ὅν... ay, / d£, 1 Α 4- 
ο SRA O Y +y 1 
Nt -x -y (Wa -W,) Vt2~x2-y tt-x* 
ΒΑ |t =v, E, (Wp-Wa) 
/ dy, ων. -f ay, / dé, ^ ps 
ο 1 O Y*y, 1 
(Ad) 


where R}? = ES - (y - y,)* and the bars on the integrals signify that 
the finite part of the integral is to be taken in the sense defined! in 
reference 10 and that the order of integration cannot, in general, be 
reversed.? For convenience set 


9v, => > En (A5) 


lFor the subsequent analysis to hold, the definition of the finite 
part given in reference 10 is essential. This definition differs from 
that given by Hadamard when it applies to multiple integrals. 


“Since the order of integration plays an important role in the fol- 
lowing development, integration first with respect to x and then with 
respect to y will be denoted fay fax f(x,y) while integration first 
with respect to y and then with respect to x will be denoted 


f dxf dy f(x,y). When the notation fff(x,y)dydx is used, the order of 
integration is immaterial. 


- ———À - -- 
A A A A αμ As RR A — c cdi a E E END c Ar 
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where In is the nth integral group on the right-hand side of equa- 
tion (Ah). 


Consider the first of these integral sets. Using equation (Al), 
we can write 


y t 
ξιᾶξ 
n= f yt 5 cro $i ax 

y-t Υ-Ψι Lē S (y-y1) ] X-a τς ιά 


= νε 2 
v-4 (x-x1)^* $i Wy (x 4Mt, y, )dt, 


| af (t-t,)* m (x-x,)* = ES 


In order to simplify this expression, the order of these integrals will 
be rearranged so the integration with respect to E, can be carried 
out first. The technique of changing the order of repeated integrals 
with strong singularities set forth in reference 10 will be used here. 
Consider the change of order in the £,, x, plane. Pretend for the 
moment, that the t, integration has 
A been carried out. Then the highest 
order singularity (since wy, is 
bounded) in the δι, x, plane has the 
order 3/2 which is weak in the sense 
that no residual occurs when the 
Ju Y-I, sequence of integration is reversed. 
The top of sketch (o) shows the area of 
of integration, so immediately 


y x44 t?-(y-y4)* 
Ii -/ ανα a dx; 
y-t x-4 t*^-(y-y4)* 


A t* - (x-x4)* 


E dë 


ry, [6 ?-(y-y,) 21777 


2 2 
-/(x- + 
tx xi) zl Va (x4 {ΜΟΙ γη} dt, 


O Ν (£-t,)2- (x-x,)?-£,7 


6M 


ee — = ee A — 
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To change order in the §&,,t, plane, consult the bottom of 
sketch (o). In this case an inherent singularity exists at the con- 
fluence of the singularity lines of the integrand; namely, where 


E, = y - y, andt, =t -J(x-x,)* + $,?. The change of order can 
therefore not be performed directly, but account must be taken of the 
existence of a residual term (see ref. 10). This residual is defined 
as the difference between the two integrals taken in different orders 
over a vanishingly small region surrounding the inherent singularity 
(the region heavily shaded in bottom of sketch (o). The residual Ry 
is then, 


A (rg*€)*- (x-x4 )? 


B HE 
= 50 uc s 
Y-Yi [£,* 7 (y-y,)*] 
t- (c ,) 8 47 ο 
PI .... πο ρα 
1 125241 1. 
t-rg-€ (t-t,)*- (x-x,)?-&* -Yg-€ 


sf ο ο (Jo 
A KHI 
Υ-γι [£,2 -(y-y,) 2 197? /(t-6,)? - (x-x4)? - &;? 


where rg? = (x - x,)* + (y - y,)*. The second integral vanishes (see 
ref. 10), and, passing to the limit e—>0 in the first integral, there 
results 


wa(xyHMt-Mro,yi) κ x [vu] 
Io 2 Tg 


R; =- £ 
i 5 


\ 


where the square brackets again mean that the retarded value is to be 
taken. Thus, the integral I, can be reduced to 


dl nn — [ο] - 
: y-t x- V/t*-(y-y,) ᾽ 


, 
A a ——Huá P cna  ——— M RM — a 


— ——— -------- — -- --------- ---- --- 
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In the same way, the integral I4 can be reduced, and 


y+t xt t^-(y-y1)* 
x 
li + lo = = zf dy, n ἄκη [wu] 


r 
y-t x= t*-(y-y,)? B 


which is recognized as Kirchhoff's formula, equation (3), with an 
acoustic plan form bounded by the circle 


(x-x4)* + (y-yi)? = tf 


The reduction of the integrals Τα, I,, I5, and Τα is quite simi- 
lar, leading to the sum 


y+t xt t*-(y-y4)? y+ Bx? 


6 

TOME [vu] 

) Ta = "A dy; 7 dx; = uf dyi 
γε... 

1 ο x t (y yi) Ο 


Ο yin t^-x* o 


l = (AT) 


Examination of the limits on these integrals shows their total area 
of integration is that shown in sketch (j). But this area corresponds 
exactly to the acoustic plan form Sa for a point in region Va! Hence, 


denoting the combination of terms in equation (AT) by p(1) we can write 
simply 


DEIN ποσῶν (48) 


Va το 
(Sa ^v. 
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It now remains to calculate the integrals Iy through Το. Desig- 


nating their total effect on the potential by p(2) ; one can readily 
show (since no inherent singularities arise in these cases) that 


-ytt xt? - (yty) T 
(2) 1 t-rı uyy, wé(x,+Mt, jy, )dt, 
g dyi ax f 


ο aldea “ο [(t-t1)? - x97] ($-6,)? - 1,7 


t-r, -y+ tox? 


AJ yy, wu (x Mt, ,y, )8t, 
o [(t-t,)? - 1,214 (t-t,)? - r,? 


d a 
O 


yy, w,(x,4Mt,,y.,)dt 
TA 1. Wye / 87 (a9) 


J -x, /M [(t-t,)* ας ro^ W(t-t,)*- pa 


where m =(x-x,)* + (y+y,)%. Now let 


t-r 
Bo a lyy, wy (x4 3Mt, , y, ) t, 


~x,/M — [(t-t,)?-r, 2] / (t-t,)?- 4? 


t-r 
[ + lyy, wy(x,+Mt,,y,)dat, 
_ δα ο κα. ey 


C(xi,yi) = (A30) 


x, >0 


- κε ο µας ——— a —— -. --- -- ------ -------- 
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In terms of this expression, equation (A9) can be written simply 


(ek Pf care (ana) 
(Sed, 


where the area (Sedy,, is illustrated in sketch (κ). 


In order to give expressions for the potential in every region of 
the wing shown in figure 1, one can show that it is only necessary to 
vary the &reas over which the double integr&tion in equations (A8) and 
(A11) are carried out. This is evident in connection with the source 


1 
portion Q'* ', for in every case 


(13). 1 [Yu] 
Q''z a NA c Ody (A12) 
Sa 


and only the acoustic plan form Sg changes with the region. In the 


case of φ(3), the part of the potential due to the existence. of the 
side edge of the wing, equation (A11) can be generalized and written 


MEE - n C(x1,y,)dx,dy, (A13) 
, Se 


where the integrands are defined in every case by equation (A10) and 
only the "reflected" acoustic plan form Se changes with the region. 
The region Se is always bounded by portions of the "reflected" circle. 


(x-x,)2 + (yty) = V^ 


and the "reflected" ellipse 
p 2 
(E κ πα) + (y+y,)% = tq 


Figure 2 shows sketches of both Sc and Sa for all regions in figure 1, 
The absence of a sketch indicates that the corresponding integral does 
not exist for that region. 
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APPENDIX B 


THE GENERALIZED INDICTAL FORCES 
The Loading Coefficient 


In order to determine total forces acting on the wing, it is first 
necessary to obtain expressions for the loading coefficient Δρ/ας. 


According to the linear theory 


Ap o9 
dS ^ Uo δε (35; 


So 1t is necessary to differentiate each of the expressions for poten- 
tial. As an example, consider, as in Appendix A, just region Va of 
figure l. The loading coefficient will be divided into two parts 


Ap /α. and Δρί5) /q,, to correspond to the potentials p(1) ana p(2), 
Thus, using equation (Α11) 


(2) -y*t xe t^ (γην): 


5) = Y J dy / oc dX4 = 
do * 9 | 
Y dal fa: xe “' 
“yt t^-x? 5 
9 X- t7- (y+y,) 
ΤΝ D O 
/ daya / i 2C ax, (B2) 
» X, (y+y,) | 


since the derivative passes the x,,y, integration without effect. 
Referring to equation (A10) for the function C(x,,y,) we next find its 
derivative with respect to t. Write T = t-t,; then for x, <0 


tex, /M λ/ lyy. W(x + Mt - MT,y,)dT 


C(xi,yi) = 
acid MEE (12-262) Jr? 2? 


τα 


----- ee -- ----- 
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and 


S e : GOES 
re Py δε porte) y at (B3) 


Ti 


Notice that if wy does not depend on (x, + Mt,) the integral term in 
equation (B3) vanishes, while if it does, then the integrated term is 
zero. Next, for x, >0, 


t y lyy, wy(x, + Mt-MT,y,) 
C(x4 ,y1) = / ———— ee T 
(12 - 162) 72 - m, 2 


n. T -rı 


and 


/ ἔργ, 9. Mt - Mr, } 
ὃς N kyy, N Ya WulX1,Y1) k ώς ot [nen 1,71) 


4 dr (mh) 
δὲ (tér?) (420.2) Vtr 8 " (72 - ro?) / 78 - y? 


In this case, both terms exist unless wy is not a function of (κ. +Mt,), 
in which case the integral vanishes. 


Substitution of equations (B3) and (Bh) into equation (B2) will now 
yield an expression for the loading coefficient corresponding to the 
influence of the side edge; 
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-y+t xt t£- (y+y, )2 
^ N byy, xy? yP dx + 


dy, 
x-Vt=-(y+y,)* (t2-ro?) y t-r,” 
-y-t x t*-( y+y,) i 


o 

E 

Lm B 
i 

oa 


Va 


ο x-4 t^-(ysy4)* 


-y* t^ -x° 
dy, 


t-ri 
a/ lyy, (x, + Mt, oy." di; 


/ [(1-ἑ1})5-αρΞ1ν (6-4, )8-r,? Yo 


ὃ ; -y4 t?-x? 
WV, X, ΥΣ ἄπ — Ml J 


ER 


— (t?-rg*) V tér, 5 ο 


dy 


t-r 
ο 1 tel 
a heyy, (x,+Mt,)"""y Pat, 


A + 
M ο σον 


-y+ t?-x? ο 
M1 [ ayı [ dx 
9 Xi (y*y4) 
t-r 


ls 
* yy, (x4 3M4t,) γα Paty 


-(x,/M) [(&-t,)9-ro? W(t-t4)*-7,* 


(85) 


The explicit form of wy, given by equation (2), has been inserted and 
it is assumed that 21. 


EE rs. Es et a at a a pa a τί 
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1. 
The portion of the loading coefficient corresponding to o ) can be 
found readily and is = 


(2 
do 


(1) 2870 [" , (εν t2-(y-y,)9) ο ο ποσο 
Va Moct? | "ο y t®-(y-y,)* 


dy, + 


1. 


tel 
y+t ἜΝ t^-(y-yi)* [x +M(t-r 
in f y By. [ l [ 1 ( ix) I dx 
r 
ο x-Vt"-(y-y,)* i: 
y+ 
M Fl yp y, --------- a, + 
; Ο 
z X- y t2- (y-y 


ΕΞ -x? l-1 
Ὥ O [x1*M(t-rg) ] 
" deg MEL 


ro 


MUT ο [x Μοτο) 


ο 
Xi (y-y4i) 


Mt? -x? 
f τα Gel er y 
e i N t2- (y-y4)* 


1 (B6) 


It is clear that, even for small values of the indices land n, 
the required integrations for the determination of total forces on the 
wing pose formidable problems. There is, however, a property of the 
loading coefficient corresponding to vertical velocity distributions of 
the type chosen here (eq. (2)) that will materially shorten the requisite 
labor. This may be expressed as follows, adopting the convention that 


Ap'™/q, corresponds to a downwash distribution proportional to 


(x4Mt)! y’: 


ο ον τον, οσο (ΒΤ) 


"M 
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ΟΥ, 


in x l-1,n 
Ap Ap 2 
a a A —— — (x,,y,t)dx,, 1-0 (58) 
Ο q 


-Mt M 
Details of Evaluating the Generalized Indicial Forces 


In calculating the generalized indicial forces by means of equa- 
tion (36), it has been shown that only the value zero need be taken for 
the index wl. Thus we must find 


ο- Μο 


on 
Fo = LE I (x+ Mt)? dx [ ye OP ay (B9) 
g Ώου 61 PS A do 


The values of the loading coefficient Ap /q. are found by differentiat- 
ing the expressions for potential given in the first part of this 
appendix. 


It is convenient, in evaluating equation (B9), to consider the inte- 
gration with respect to y first. Setting 


8 & on 
EJEA (B10) 
ο i “o 


it is found that L seems to have different representations according 
to the interval in which x lies. These expressions can, however, all 
be η the same formula. The portions of L corresponding to 


the parts 9 and q 2! of the potential are similarly signified, and 
we have 


DL. . Fon ἯΝ 1 (-19 E όν | Ko(mte) + Kyle) | - 


πομοΐ 8 (n+g+1)! 
[n/2] (5) Pere 

2 ) (a Lu [ren + (0-2) | } (B: 
ΗΞΟ 


SHE | Klute) + (ae) | | (B: 
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οος-1(-κ/9) 
Ko(ntg) = t3 ΕΥ" R.P. n sinb*B*i gag 


cos”1 (Xy /tm) ees 


6de 


Ky(n*g) 


J(n,g) 


LO. Un j: ae τ 


y 
a [ro 


[T 


and [n/2] means the greatest integer contained in n/2. The function 
J(n,g) may be expressed as summations, and it has the property 


J(n,g) = J(g,n) (B13) 


The sum formula is, with g+ p= 


[p/2] 


E 28+1 p-2i+1 2g+2 
en y E A e]. 
y 8-1 T g-1 
(-1 ER j {2113 1 p-21+23+3 28-23-1 
a 21) ) CU (Ea, E j C mr zn 
izo p 


y RE 


), (8) ye ES 23+1 , 83) ES 2 | (814) 
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Values of the function J(g,n) 


2329 5 - 
1638} 


Tt + 


"65536 693 


where $4 18 the binomial coefficient 


u — c - 
ei (21)! (p -21)! 
and Blp,q) is the beta function 


7 -=L 
B(p,a) -f 22 1-17 ax 


O 


0/2 2a-1 
= El gin PL € cos?” gae (B15&) 


O 


= Τ(ο) r (q) /T(p + a) 


a ee a A A mm mr TR 
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The function J(g,n) has been calculated for g,n taken 0,1,2,3,4,5. 


Because of the property (B13), it is only necessary to give a triangular 
array, which appears in the above table. 


Now consider the functions Κο(ν) and Κμίν), defined after equa- 
tion (B12). It is convenient, for computational purposes, to express 
these in terms of the incomplete beta functions, defined as 


-1 
cos *(x) 
B, x2 (ρ,α) = 2 n sin’? lg cos + - edo 
O 


(815b) 


1-x* 


/ gP-2 η 69171 at 


ο 


A tabulation of the incomplete beta functions is available in refer- 
ence ll. Note that when the symbol B is written without a subscript, 
the complete integral is meant, that is, in equation (Bl5b), x equals O. 
It is necessary to exercise some care when interpreting Κο(ν) and Κμίν) 
as beta functions because of the upper limit. Thus, since 


-1 
οος-τ(-κ/ 9) 
Κο(ν) = ο RP. f. sin '* eae 


O 


we have the following cases: 


(i) x >t, R.P. cos”1 (- x) =x 


Κο(ν) = 4VtL (=, 2) 


E x 


κοιν) = £57 E C i) πο 3) | 
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(iii) -t € x € 0, R.P. cos”? (- E = 


23 


cos 1(-X 
i t 


Ko(v) = “η AC 2) 


(iv) -Mt <x <- t; R.P. cos! (- E) 
Ko(v) = 0 
A similar line taken with Ky(v) leads 
(1) x >t, Ky(v) =0 


_i SLM, νει 
(1i) M <x t, Ky(v) =5 pm. |» T 


to 


anfen) Ἔ» B | 


Sy -(25p/tm) Es, =) | 


(iv) -Mt <x € -t, Ky(v) = tnt B( 42 5) 


ο 7 9 


The generalized indicial force pon can now be expressed as 


JE 


a J(g,n) 

ο... ρα . , n 

F = ol- 
jg MU citen uus τς | ατα parm 7 (-1) 


[n/a] 


xcd ,8tn*i-ep 
n em]- 5 (δι)  — 
M ΕΗ} ginil-2n 


μ-ο 


ο 2 
a te me Pg e Pao ——— - —— 


e k (g+n) + 


k (24-1) + ας (eu-1) |) (25) 


—A-—— ------- ee — — —— — - 
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54 
where 
: c-Mt 3 cos “(-x/t) 
"Io (v) = + (x4 Με)". dx Lis R.P. J sin”? cao | (B17) 
"Mt A 
c-Mt cos  (xm/ πι) 
ges (v) =f (x4 Mt)? αχ] Ἡ tm s RP. f sin t+ 090 | 
-Mt ο 


(818) 


It is convenient to express these forces in terms of dimensionless 
quantities. Thus setting 


«cb owl 
we have 
1-Mto 
+V+2 
κτὸ (y) - οὗ (xo 4 Μορ) dxo ran R.P. 
-Mto 
cos71(-xg/to) | 
l et ono |= VE τοῦ (y (519) 
O 
3 ἣν πα. V+1 
x 4442 +Mt 
I (Y) = c r (xo + Mto) J axo (ο) Β.Ρ. 
-Mto 
+ ν TVA 
Sotto sip MH gag | = c? Ty? (v) (820) 
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and 
J(g,n) co 
S| Ἕθηα- ogtn 9 
(-1) 2 ME | τὸ (gin) + -of 
? 
(n+g+1)! mi, 
[n/2] A ) (i) O< b< IMI) 
j ( l n ο 
gin) |. bis CNN 
M Εμ) σιαεῖ-ρμ 
H=0 lo 
| τ (20-1) + τὸ (en-2) | } (B21) 
ο IM 
The integrals τὸ (ν) and τὸ(ν) M? 
can be simplified by reversing the /- Mio 
order of integration. This can be -h 


accomplished in a straight-forward 

manner by merely inspecting the region 

of integration in the X9,0 plane. 

Consider first the integral τ (v). |. -Mt 
Depending upon the relation between the 7 E 
chord length and the time, we see - (i) WMH) < 6 <//M-1) 
from sketch (p) - that πμ the 
order of integration results in three 
different possibilities for the upper 
limit of the € integral. However, if 
we define X, such that 


Xo 


(i) Xa = to; O< ος t. 
L L 
( ) Xo = ο. Mel 9 Cl A 
/-Mfo 
1 
(iii) X. = -t,; —— et 
O ο M-1 = O 


-Mth 


(iii) IAM-l) fh, 
Sketch (p) 
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then, in every case, τὸ (v) can be written 


--(- 4-1 | 
τον ΤΙ cos" (-Xo/to) 4 ote i 
τὸ (ν) = e / sin"'* gag - "y κ} y 
ο 


I-O 


sin '* O cosTedo (B22) 


re 


O 


and, similarly, it can be shown that 


τ. GQ) οσο τν Fi O0/"O sn * 6d0 + 
- 115 


gonia οος” τ(-Χο/ το) 


μα” » (-1) n et sin". 8 cos” 0d0 
(B23) 


BM 


‘ated with the subscript 1 have 
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APPENDIX C 
DERIVATION OF RECIPROCITY RELATIONS 


According to reference >, the reciprocity relation for general 
three-dimensional unsteady motion can be written 


Ap, 
"i ᾳ (xi,yi,t)Wa(xi,y,,t,) dx,dy,dt, = 
ο 
V 


Ap 
ni 7m (X59 52 t2) W1 15,3556.) dx,dydt,, (C1) 
ο 
V 


where the volume of integration V is that swept out in x,y,t space by 
the wing. The subscript 1 refers to the wing moving in the forward direc- 
tion and subscript 2 refers to the wing moving in the opposite direction 
in the same manner. The coordinate systems are related by 


Xy = -Xptc-MT 

Ji" “yat 28 

Ὅ =- to+ T 
where s,c are wing semispan 
and chord, respectively, and T 
is some fixed value of time. 


These quantities are elucidated 
in sketch (q). 


Now let the wing associ- 


the vertical velocity distribu- Sketch (q) 


tion 
X, + Mt, t S - yi B 
"Gy t) (2) (= 


and that associated with the subscript 2 have 


Xo + Mt, J /s -y> NE 
WolXo)Yp1to) = (18.339) (== 


nn μμ. A A ον g A e «ας A 


r— ————— —————-———A———— 
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7 Xo+ Mt, N {7ο - sn 
(1 i tris y (a= ¥ 
e ο 


Substitution of these results into equation (C1) yields 


Then 


Wo (X1 31581) 


m -Mt 
S i X4 + MG, J E" Jz - 8 6 Ap?” 
dt, dX4 1 - BF n dys TA qo E 
O -Mt O 
Y cone Mt Jg 
Xo - 
{ας f mx, (i. TETTE ϱ) fm (ες E- 
ο 


Μορ (C2) 


Equation (C2) can be differentiated with respect to T, yielding 


MT 
ux FECE 
ο 
c-MT 7 2 n j 
MT yo -8 Apv& 
{ πα πο "ο T 
ο 


-MT 


The binomial expansion is now performed: 
in 


ο Ε- ul Fa (an) ` 
aaa E ὦ 


Ho 


5 


— — 
See a tr ee παν ο ag gt gpg et 
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In equation (C3) the spanwise integration is carried over the whole wing, 
but it can easily be reduced to integration over, say, the left panel by 


use of the factor [1+ (-1)8" 9/5. Thus, equation (C3) can be written 


cas cap (2) Bette ES 


ΗΞΟ -MT 
7 
8 | E 1 18407 Jo 
PAR a Y NAAA 
do se 
O j μξο 
c-MT p 8 x 
f e) f ev (T8 Apd8 
ο ο 
-MT O do 


By comparison with equations (36) and (37), it is seen that the integral 


terms in the last equation correspond to the generalized indicial forces 


i 
f vs and "S So that the summations can be written 


γω’) E y C (1) eie (ch) 
H=0 ΗΟ 


where the quantity (g+n) must be an even number. 
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Figure l.- Regions used in the analysis of a rectangular wing in supersonic unsteady motion. 
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Figure 2.- Sketches of areas of integration, Se and Sa, for all regions 
in figure l. 
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Figure 2.- Continued. 
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Figure 2.- Concluded. 
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